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INTERPRETING POWER ANISOTROPY MEASUREMENTS IN PLASMA TURBULENCE 
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ABSTRACT 

A relationship is derived between power anisotropy and wavevector anisotropy in turbulent fluctuations. 
This can be used to interpret plasma turbulence measurements, for example, in the solar wind. If fluctuations 
are spatially anisotropic, then the ion gyroscale break point in measured spectra in the directions parallel and 
perpendicular to the magnetic field would not occur at the same frequency, and similarly for the electron 
gyroscale break point. This is an important consideration when interpreting solar wind measurements in terms 
of anisotropic turbulence theories. Model magnetic field power spectra are presented assuming a cascade of 
critically balanced Alfven waves in the inertial range and kinetic Alfven waves in the dissipation range. The 
variation of power anisotropy with scale is compared to existing solar wind measurements, and the similarities 
and differences are discussed. 

Subject headings: magnetic fields - magnetohydrodynamics (MHD) - plasmas - solar wind - turbulence 



1. INTRODUCTION 

Plasma turbulence is observed to be anisotropic with re- 
spect to the magnetic field direction. For example, in the solar 
wind, the observed power and scaling of turbulent fluctuations 
vary depending on the an gle between the local mea n field 
and the sampling dire ction (iBieber et al. 1996; Horb ury et alj 
l2008t iPodestal 120091: |Osman & Horbury 2009). Correla- 
tion functions have also been observed to b e anisotropic in 
the solar wind (Crooker et al. 1982; Matth aeus et al.l 119901; 
lOsman & Horbu rv 2007; Wevgand et al. 2009) and laboratory 
measu rements (Robinson & Rusbridge 1971; Zweben et al. 
[19791) . 

Recent theories of plasma turbulence assume anisotropy 
jGoldreich & Sridhar 1995; Boldyrev 2006; Galtier 2006; 
Lithwick et al. 2007; Gogoberidze 2007; Chandran 2008i 
Beresnyak & Lazarian 2008; Podesta & Bhattachariee 2001 
Schekochihin et al.l 12009 ) and anisotropic energy trans- 



fer has been seen in simulations ( S heb alin et al .1 1198 31; 
Cho & Vishniac 2000; Maron & Goldr ej^hn2001l: ICho et al] 
2002; ,Cho & Lazariani,2004. ,2009) . The theories usually de- 



scribe the anisotropy in terms of the fluctuation wavenumbers 
parallel and perpendicular to the mean magnetic field direc- 
tion, and k^. For example, Goldreich & Sridhar (199^ 
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used the "critical balance" assumption to obtain k\\ ^ for 
magnetohydrodynamic (MHD) turbulence and, more gener- 
ally, theories often assume k±^ ^ k\\. In the solar wind, how- 
ever, it is the anisotropy in power at a fixed scale that is 
often measured for practical reasons, rather than the spatial 
anisotropy of the fluctuations. 

In Section|2] the relationship between power anisotropy and 
wavevector anisotropy is derived. A critically balanced model 
is then presented in Section |3] to illustrate that a break point 
in an anisotropic spectrum may occur at different scales when 
the reduced spectrum is observed in different directions. The 
implications of these two sections for recent solar wind mea- 
surements are discussed in Section|4] 



2. POWER ANISOTROPY AND WAVEVECTOR 
ANISOTROPY 

The correlation function of a turbulent field, for ex- 
ample, the magnetic field, B, can be defined as C(x) = 
(B(r-l-x) • B(r)), where the angular brackets denote an ensem- 
ble average over positions r. The three-dimensional energy 
spectrum can then be defined as the Fourier transform of the 
correlation function, E(\i.) = J C{k)e~'^'^d^x. 

A single spacecraft in the solar wind measures the turbu- 
lent field as a function of time, B(f). Since the solar wind 
velocity, Vsw, is much larger than the wave speed (often taken 
as t he Alfven spe ed in the inertial range), Taylor's hypothe- 
sis (lTaylor|[T938[) is usually well satisfied, meaning that the 
measured time variations correspond to spatial fluctuations 
in the plasma. Ax = -VswAf. Because a single spacecraft 
gives a one-dimensional cut through the plasma, the full three- 
dimensional spectrum cannot be measured but instead a re- 
duced version is obtained (Fredricks & Coroniti..l976.) . This 
reduced spectrum, defined as P(k) = J E(k')S(k-k' ■ Vsw)'^^k', 
where Vsw is the solar wind direction unit vector, is the three- 
dimensional spectrum integrated over the directions perpen- 
dicular to the measuring direction. Assuming axisymmetry 
about the magnetic field, this reduced power spectrum also 
depends on the angle, 6, of the field to the one-dimensional 
measurement direction, P(k, 9). In Cartesian coordinates, this 
can be written as a dependence on the parallel and perpendic- 
ular wavenumbers, P{k^^,k±). 

Figure [U is a schematic of reduced power contours with re- 
spect to the parallel and perpendicular wavenumbers. Power 
anisotropy is usually measured at a fixed scale, indicated by 
the red dashed line which is at a fixed radius from the ori- 
gin. At different points along this line, a different reduced 
power is sampled; this effect is readi ly seen in the solar 
wind (Bieber et al.''1996'; 'HorburiLeLalJ|2008'; Podesta 2001 
Osman & Horbury 2009). A relationship between power 
anisotropy and wavevector anisotropy will now be derived. 
Note that the derivation does not depend on any particular 
contour shape; the elliptical shapes in the figure are for illus- 
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Parallel Wavenumber, Ij^ 



Figure 1. Schematic of reduced power contours (solid blue lines) elongated 
in the field parallel direction. Power anisotropy measurements are made at a 
fixed scale (dashed red fine). 



trative purposes only. 

Let us consider two contours of size and 
ik^\2,k±2) such that k\\2 = (Figure[T]i- We will assume that 
power anisotropy is being measured at this wavenumber, k = 
k[\2 = k±i. Let reduced power in the parallel and perpendic- 
ular directions be defined as P|| = P(k\\,0) and P± = P(0,k±), 
and a be the scaling exponent in the perpendicular direction, 
P± ~ k'^. Dividing the equations for P± for each contour we 
get 



P±2 



k±2 
k±i 



(1) 



By the definition of a contour, P±2 = P\\2 ™d from the above 
definition of the two contours, kp = k±i. Substituting these 
into Equation ([T]i gives 



P\\2 



C||2 



(2) 



Since P±i/Pp is just the power anisotropy at fixed wavenum- 
ber k (as in solar wind measurements), the numeric subscripts 
may be dropped. Rearranging Equation (|2]i we get 



(3) 



This relationship is independent of the scaling of the parallel 
spectrum and allows us to calculate the wavevector anisotropy 



of contour 2 from a measurement of P 
tionship can be derived for contour 1 



'±/P\\- A similar rela- 



(4) 



where f3 is the scaling exponent of the parallel reduced spec- 



trum. Pi I 



-|| ~ /C|| . Although it is possible to infer the wavevec- 
tor anisotropy from the interpolation of measurements such as 
Figure 1 of Horburv et al. (2008), the relationship given here 
allows it to be found from the power anisotropy, a quantity 
more easily measurable in the solar wind. 

The k± and ^|| of turbulence theories usually describe typ- 
ical length scales associated with the fluctuations. It is usu- 



ally assumed that second-order stati stics, such as power, re- 
late to these quantities, for example ICho & Vishniad (l2000h 
state that in their simulations, contours of second-order struc- 
ture functions "reflect the shapes of the eddies." Under this 
assumption, the wavevector anisotropy of power contours can 
be thought to describe "typical" wavevector anisotropy of the 
fluctuations. 

3. FORM OF THE CRITICAL BALANCE REDUCED 
POWER SPECTRUM 

The "critical balance" assumption states that in a turbulent 
Alfven wave (AW) cascade, the linear wave timescale and the 
nonlinear energy transfer timescale are comparable. It was 
introduced explicitly by Goldreich & Sridhar (1995) and an- 
ticipated in the work of Higdon (1984). When applied to 
inertial range MHD turbulence, the spectral index of the re- 
duced spectrum in the perpendicular direction is -5/3 and the 
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wavevector scaling is k\\ ^ kj_ . A reduced spectral index of 
-2 in the parallel direction follows from these statements. 

There is evidence in the solar wind inertial range 
for both the Alfvenic nature of the turbulence (e.g . 
iBelcher & Davis|[T97Tt iHorburv et al.l[T99l iBale et allllOOl 
and the anisotropic scaling dHorburv et al.l |20()8Ik It ap- 
pears that this scaling is only detectable when observ- 
ing with respect to the scale-dependent lo cal mean mag- 
netic field (e.g. Cho & Vishniac 2000; H orburv et al.1 120081: 
iBeresnvak & Laz arian 200^, i.e., the mean field at the scale 
of each fluctuation being measured, rather than a global 
large-scale average field. Although MHD is a fluid theory, 
ISchekochihin et al.l (|2009) have shown that reduced MHD, 
an anisotropic limit of MHD containing the Alfvenic fluctua- 
tions, can be derived for a collisionless plasma at scales larger 
than the ion gyroradius. This may explain why the MHD scal- 
ings are seen in the collisionless solar wind. 

At scales smaller than the ion gyroradius, commonly 
termed the "dissipation ran ge," there is evidence for ki- 
netic Alfven waves (KAWs) (|Bale et al.l l2005': 'Sahraou i etaP 
|2009l). These are linear modes of electron reduced MHD, an 
anisotropic theory derived for collisionle ss plasmas at scales 
between the electron and ion gyroradii (ISchekochihin et al] 
120091) . When the critical balance assumption is applied to a 
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KAW cascade, the wavevector scaling becomes k\\ k^ and 
the predicted spectral indices for the magnetic field are -7/3 
in the perpendicular direction and -5 in the paral lel direction 
dCho & Lazarianll2004t ISchekochihin etani2009h . 

In this theoretical framework, the break between the in- 
ertial range and the dissipation range is predicted to be at 
k±Pi ^ 1, where p,- is the ion gyroradius, but if the fluctua- 
tions here are anisotropic then their parallel length should be 
larger, < 1. This would imply that the observed break 
points in solar wind measurements are at different spacecraft 
frequencies for the reduced spectra in the parallel and perpen- 
dicular directions. A similar effect would be expected at the 
electron break scale, k±pe ^ 1, where the difference in break 
frequency between the spectra in the parallel and perpendicu- 
lar directions may be even greater if the anisotropy continues 
to increase throughout the dissipation range. 

Since the observed break points would be at different scales 
for spectra in the parallel and perpendicular directions, a 
schematic of the reduced spectra can be divided into five 
ranges (Figure |2]i. In range 1 both spectra display AW scal- 
ing; in range 2 the spectrum in the parallel direction has KAW 
scaHng and the spectrum in the perpendicular direction has 
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Figure 2. Schematic of magnetic field reduced power spectra in tlie parallel 
and perpendicular directions for critically balanced AW (solid black lines) 
and KAW (dashed blue lines) turbulence. 



Table 1 

Power Scaling Exponent Predictions for Ciitically 
Balanced AW and KAW Turbulence 



Range 


Pj^ Scahng 


P|| Scaling 


Px//"!! Scaling 


1 


-5/3 


-2 


1/3 


2 


-5/3 


-5 


10/3 


3 


-7/3 


-5 


8/3 


4 


-7/3 






5 









AW scaling; in range 3 both spectra display KAW scaling; in 
range 4 the spectrum in the parallel direction is below the elec- 
tron break scale and the spectrum in the perpendicular direc- 
tion has KAW scaling; in range 5 both spectra are below the 
electron break scale. Predictions for fluctuations smaller than 
the electron scale do exist but have not been included here. 
Gyrokinetic theory predicts scalin gs for an e lectron-entropy 
cascade, valid for k±pe^ I (.Schekochihin et al. 2009 ). how- 
ever it has been suggeste d that it is no t applicable to the solar 
wind in this range ( Howes et al.ll2008l) . In Figure El the (log- 
arithm of the) power anisotropy can be thought of as the ver- 
tical distance between the spectra and the (logarithm of the) 
wavevector anisotropy as the horizontal distance. 

The scalings for each of the ranges in Figure|2]are given in 
Table [T] Also listed is the scaling of P±/P\\, which follows 
directly from that of P± and P|| and is shown in Figure [3] 

In the inertial range P±/P\\ scales as A;'/^ which steepens to 

^10/3 when P|| reaches the ion break scale and then becomes 

shallower at k^^^ when P± reaches the ion break scale. 

The width of the KAW range in the P± spectrum is pre- 
dicted to be pi/pe = ^/Tiny/^]me, where m,- and nig are the ion 
and electron masses, and for the model spectra in Figures |2] 
and[3]the temperatures have been assumed equal, T] = T^. Al- 
though Ti/Te is o f order unity in th e solar wind, there is some 
variation ( Brun o & Carbondl2005h so the extent of the possi- 
ble KAW range may vary. The width of each of the ranges 
in Figures |2] and [3] also depends on the amount of anisotropy 
present. For example, as the anisotropy at the ion break scale 
increases, the size of range 2 increases but the size of range 
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Figure 3. Schematic of reduced power anisotropy as a function of scale for 
critically balanced AW and KAW turbulence. 

3 decreases. It may even be the case that if the anisotropy 
is very strong at the ion break scale then range 3 may not be 
present at all, meaning it would not be possible to measure 
KAW scaling in and P|| at the same frequency. 

One of the assumptions used when constructing these 
model spectra was that they c ontain no additional ener gy in- 
jection or loss. As noted in I Schekochihin et akl (12009 ). at 
the ion break scale some energy may be transferred from the 
Alfvenic cascade channel to a purely electrostatic "entropy 
cascade." The amount of energy transferred, if any, is un- 
known so cannot be included in our model spectra here. It 
is also possible that power may be injected into the cascades 
from other sources such as plasma instabilities, for example, 
the firehose and mirror instabilities, which evidence suggests 
may be important in the solar wind (Bale et al. 2009). Includ- 
ing effects such as these in this model is beyond the scope of 
this Letter 

It should also be noted that only the Alfvenic part of the cas- 
cade in the inertial range is dealt with here. This is relevant 
to the solar wind, which is primarily Alfvenic in nature (e.g. 
'Belch er &DavisllT97Tl; iBale et all BOOl iBruno & Carbong 
2005) and in which any compressive (non-AIfvenic) fluc- 
tuations are, on theoretical groun ds, not thought to in- 
terfere with the Alfven ic cascade (ICho & LazarianI 120031 : 
ISchekochihinetani2009l) . 

4. COMPARISON TO SOLAR WIND MEASUREMENTS 

The only published measurement of the variation of power 
anisotropy, in the solar wind across both inertial and 

dissipation ranges, which we are a ware of, is in the lower 
panel of Figure 7 of Podesta" (2009). The scaling arguments 
in Section [2] can be applied to these measurements to obtain 
estimates of the wavevector anisotropy at various scales. For 
example, at the high-frequency end of the inertial range, at 
« 0.2 Hz, P^lP\\^ « 7 which, using a = 5/3 in Equation (O, 
means k^lk\\ « 3. A value of a = 5/3 was used for this cal- 
culation since it is the prediction of critical balance MHD tur- 
bulence (Goldreich & Sridhar 1995) and i s close to the mea- 
sured value of 1.65 for this data interval ( lPodestall2009 ). In 
general, break points in spectra seem to have a rollover rather 
than a clean break in scaling so this result is approximate. It 
is also possible that since P\^ is measured in the bin 0-6° and 
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not at exactly 0°, values of P±/P\\ and therefore k±/k\^ may 
be underestimates. The value obtained here, kj_/k^\ « 3, was 
used to set the anisotropy at the ion break scale in Figures |2] 
and[3] 

Figure [3] of this Letter can be compared to Figure 7 of 
IPodesta (.2009) in which the frequency corresponding to kp, = 
1 is « 0.5 Hz. In the inertial range, both graphs have a shal- 
low slope with the measurement steeper than the prediction, 
although the uncertainties in the measured slope may be sig- 
nificant as discussed in Podesta (2009). Between 0.4 Hz and 
1 .0 Hz, the measured power anisotropy increases with a steep 
slope. The scales at which this happens approximately corre- 
spond to ranges 2 and 3 in Figure|3] It may be the case, there- 
fore, that the steep slope of P±/P\\ seen in dissipation range 
solar wind measurements is due to critical balance scaling, in 
particular the steep scaling of the KAW P|| spectrum. 

One difference between the figures is that the decrease in 
power anisotropy between 0.2 Hz and 0.4 Hz in the r neasure- 
ment is not present in Figure |3] As discussed in iPodestal 
(120091) . this is caused by an increase in the parallel power 
and may be due to parallel waves, for example, from plasma 
instabilities. For frequencies above 1.0 Hz the measured 
anisotropy decreases, which is due to the flattening of the f || 
spectrum. The scale at which this begins (1.0 Hz) is close 
to the predicted electron break scale for the P|| spectrum, al- 
though without knowing the electron gyroradius, the exact lo- 
cation of this is not clear One possibility, therefore, is that the 
cause of this decrease in P±/P\\ may be due to the P|| spec- 
trum flattening above the electron break scale. Another possi- 
bility is that the cyclotron resonance may have been rea ched 
here, causing a change in behavior dHowes et al.ll2008l) . At 
these high frequencies, however, measurement effects, such 
as magnetometer noise, may be important and one must be 
cautious when drawing any conclusions from this range. 

Extrapolating the model spectra to larger scales, it can 
be seen from Figures |2] and Q\ that power (and wavevector) 
isotropy is reached at kpi w 10"^. This corresponds to a 
length around 10'' km or an observed spacecraft frequency 
of 5 X 10"^ Hz. This is close to the observed break be- 
twee n the low-frequency power law and the inertial range 
(e.g. lBavassano et aLlll982l:lBruno & Carbone 20 051) and also 
the solar wind correlation length (e.g. Matthaeus & GoldsteinI 
[1982), scales usually associated with the outer scale of the 
turbulence. 

The solar wind sp ectrum has recently been observed at 
scales near kp^ = 1. lAlexandrova et al.) (120091) suggest that 
just above this scale there is an e xponential falloff in the spec- 
trum, and ISahraoui et al.l (120091) suggest that below it, there 
is a further steeper power law. Both of these studies involved 
solar wind intervals where the magnetic field was not aligned 
with the solar wind direction so that one woul d expect to see 
the spectrum in the perpendicular direction. ISahraoui et al.l 
((2009) also plot the spectrum of the parallel component of the 
magnetic field, which, it should be pointed out for clarity, is 
not the same as the spectrum in the parallel direction. 

5. SUMMARY AND CONCLUSIONS 

In Section |2] a relationship was derived that allows the 
turbulent wavevector anisotropy, k±/k^\, to be inferred from 
power anisotropy measurements, P^/P\\. This is independent 
of any particular turbulence theory and only assumes power- 
law scaling in the parallel and perpendicular directions. Using 
this relation and existing solar wind measurements (Podesta 



12009!) the wavevector anisotropy near the ion break scale was 
estimated to be A:^ « 3, although this may be an underesti- 
mate due to the finite angular resolution of the measurements. 

Model spectra of critically balanced AWs (for the inertial 
range) and KAWs (for the dissipation range) were presented 
to illustrate that break points do not occur at the same scale in 
the observed spectra in the parallel and perpendicular direc- 
tions if the turbulence is anisotropic. The variation of P±/P\\ 
with scale was calculated from these model spectra resulting 
in five ranges, three of which have predictions of how P±/P\\ 
scales: a 1/3 range, a 10/3 range, and an 8/3 range. If the 
wavevector anisotropy is significant then some of these ranges 
are small and may not even be pres ent. Some of the se features 
can be seen in the measurements of IPodestal (l2009h . The main 
difference is the extra parallel power at the ion break scale 
seen in the measurements, which may be due to energy injec- 
tion mechanisms at the ion gyroscale. 

Although critically balanced AWs and KAWs were used 
in Section |3] the ideas also apply to anisotropic theo- 
ries of plasma turbulence in general. Some of these, 
for e xample, are a nonlocal cascade model (GogoberidzS 
'2007), turbulence with dynamic alignment (Bo ldvrevI 
2006; Podesta & Bhattacharjee 2009), wave turbulence in 
Hall MHD (jG altier 20 06|), and imbalanced turbulence 
dLidiwick et alT2007; Beresnvak & Lazarianl2008l: IChandranI 
120081) . The anisotropy relationship derived in Section |2] and 
the observational considerations discussed in Sections [3] and 
|4]are also applicable to these theories and any possible exten- 
sions of them into the dissipation range. 
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